EPJ manuscript No. 

(will be inserted by the editor) 



o 
o 

(N 
-i— > 
o 

O 

o 

m 



Steps and facets at the surface of soft crystals 



P. Nozieres 1 , F. Pistolesi 2 , and S. Balibar 3 



1 Laboratoire d'Etude des Proprietes Electroniques des Solides CNRS, B.P. 166, 38042 Grenoble Cedex 9, France 

2 European Synchrotron Radiation Facility, B.P. 220, 38043 Grenoble Cedex 9, France 

3 Laboratoire de Physique Statistique de l'Ecole Normale Superieure associe aux Universites Paris 6 et Paris 7 et au CNRS 
24 rue Lhomond, 75231 Paris Cedex 05, France 



Received: 26 June 2001 



> 

CO 

a 

C 

o 
o 



Abstract. We consider the shape of crystals which are soft in the sense that their elastic modulus fi is small 
compared to their surface tension 7, more precisely \xa < 7 where a is the lattice spacing. We show that 
their surface steps penetrate inside the crystal as edge dislocations. As a consequence, these steps are broad 
with a small energy which we calculate. We also calculate the elastic interaction between steps a distance 
d apart, which is a 1/d 2 repulsion. We finally calculate the roughening temperatures of successive facets 
in order to compare with the remarkable shapes of lyotropic crystals recently observed by P. Pieranski et 
al^^. Good agreement is found. 



PACS. 81.10.Aj Crystal morphology and orientation - 
Liquid crystals 



61.72.-y Defects and impurities in crystals - 61.30.- 
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1 Introduction 

An clastic material is characterized by an elastic modu- 
lus /z which is an energy per unit volume. Any interface, 
whether with a vapour or a liquid phase, has a surface 
energy 7 per unit area. The ratio 7//Z is thus a charac- 
teristic length which should be compared to other lengths 
in the problem, for instance the lattice spacing a. More 
specifically, let us define 



d = 



(l-<r)7 
A* 



(1) 



where a is Poisson's modulus. Three situations may occur: 
(i) For usual materials in contact with their vapour, do is 
smaller than a. Indeed, according to Frank and Stroh|J, a 
usual value for 7/// is 0.1a and Poisson's modulus is usu- 
ally equal to 1/3. As a consequence, we will see that steps 
are edge dislocations bound to the surface. The familiar 
"terrace, ledge, kink" ("TLK") modelf§§ holds, (ii) At a 
solid-melt interface 7 is smaller because the solid density 
is not very different from that of the liquid. Then do <C a: 
elastic strains due to surface pertubations like steps are 
only small corrections. 

(iii) It may happen that d Q is larger than a : then it pays 
to spend elastic energy in order to gain surface energy. As 
we shall see below, steps are dislocations which are buried 
a depth h below the surface. A new length scale appears 
in the problem. The steps are broad and their energy is 
strongly reduced by an elastic relaxation. That should al- 
ways occur if the material is soft enough. A spectacular 



illustration of this situation is given by the lyotropic crys- 
tals recently studied by Pieranski et al. These crys- 
tals are water solutions of non-ionic surfactant molecules 
which arrange in a loose cubic structure with a of order 
50 A. The structure is soft (/i = 10 7 erg/cm 3 ), but the 
surface tension is only a little less than that of pure wa- 
ter (7 = 25 dyn/cm). As a result, the ratio do /a is much 
larger than in usual crystals. This is a typical example of 
the anomalous case of "soft crystals" which we consider 
in this article H . 

Pieranski et al.jjH) have observed that the facetting 
of their lyotropic crystals is remarkable: up to 60 different 
facets are present on the equilibrium shape. Our first mo- 
tivation was to understand the origin of the large number 
of facets on these crystals. As we shall see, it is the large 
size of the lattice spacing a. Pieranski et al. also invoked 
the large value of the step-step repulsion, but, as we shall 
see, it is compensated by the small value of the step energy 
in the calculation of the roughening temperatures. In sec- 
tion II, we calculate the unusual structure and the energy 
of steps in such soft crystals. In section III, we calculate 
the interaction energy of two steps a distance d apart, and 
we show that it is a 1/d 2 repulsion whose origin is surface 
capillary forces within the step profile. In Section ^ we 
calculate the net macroscopic surface energy of a vicinal 
surface, with regularly spaced steps, which is the starting 
point for a discussion of shapes. Section || is devoted to the 
roughening transition of such vicinal surfaces. In the last 
section, we use our results to calculate how many facets 
should be observed on the crystals of Pieranski et al. UP 
and we compare with their experimental results. 
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2 Structure and energy of a single step 



displacement at the surface is 



For simplicity we consider a simple cubic crystal with lat- 
tice spacing a. If no relaxation occurred, a single step on 
a (100) facet would expose an extra surface a (per unit 
length of the step). In this case, the step energy j3 would 
be 7a. Such a step is a remnant of an edge dislocation 
along the y axis with a Burgers vector along the z axis as 
it moves towards the surface. One way to reduce its en- 
ergy is to move it inside the crystal as an edge dislocation 
which is now a depth h below the free surface z = 0. A 
similar phenomenon was considered by Lejcek et al.0 in 
the different context of smectic A liquid crystals. The dis- 
location introduces an extra lattice plane on half the area 
and produces a rounded step at the surface. Burying the 
dislocation rounds off the step and decreases the surface 
energy. The price to pay is the dislocation elastic energy. 
As we shall see, this price is low for soft crystals, but it 
would be too high for usual crystals. 

In an infinite medium, the displacement uf^ and stress 

afj created by the edge dislocation are known || , as well 
as the total elastic energy per unit length 







(») 



4tt(1 



■In 



(2) 



where r Q is an appropriately defined core radius of the 
order of the lattice spacing a. i? max is the size of the sys- 
tem which will be replaced by h in the presence of a free 
surface (images suppress divergences at infinity). Qualita- 
tively we can guess the final result easily: the step width 
at the surface is of order h and the profile lengthening 
is thus w a 2 /h. Balancing f3(°> against the capillary en- 
ergy ja 2 /h yields at once h 7(1 — a) / fi = d : steps 
are buried in soft crystals! But a quantitative discussion 
requires more effort! In order to calculate this elasticity 
problem, we shall proceed in four successive steps. As a 
zeroth order approximation, we will start with an infinite 
medium, and our notations for stresses and displacements 
will have a superscript (o). Suppressing the upper half 
space modifies the strain field in the lower part, thereby 
adding an extra displacement u' 1 ' at the surface z = 0. 
The energy stored below is modified, and we must sub- 
tract the energy originally stored above. In a final step we 
must add a correction u cap due to the existence of surface 
tension. 

Since we first ignore the effect of surface tension, the 
normal stresses a z °z and have to be cancelled at the 
crystal surface. If we dealt with screw dislocations, we 
would treat the elasticity of a half space by introducing 
image dislocations. However, this method does not work 
with edge dislocations: adding a symmetric edge dislo- 
cation with opposite Burgers vector does cancel out the 
stress a zz , but it doubles cr zx ! We are thus forced to solve 
the elastic problem directly. In an infinite medium, the 
upper half exerts a force = 0^°) on the lower half. 
Cancelling this force is the same as applying an opposite 
force —F^°\ resulting in an extra displacement u[ and 
stress on top of the original uf\ 0$ . This extra 



u { p (x) 



Xij (x - x') {-Fj (a:')) dx' (3) 



where Xij ( x ~ x ') is the two dimensional surface response 
function, given for instance in Landau-Lifschitz j^] (see eq. 
A6 in the appendix) . The net normal elastic displacement 
that controls the surface profile is 



(4) 



In order to calculate the clastic energy, two corrections are 
needed 

(?) Add the energy lowering upon relaxation of the forces 



F- 



Co) 



5(3^ 



+ OO 



dx I u^dXa^(l-X) 
Jo 

+ 00 

(°) 



dxdx' Xij {x-x')a^{x)a { ;){x') (5) 



(ii) Subtract the part of 0^ that was stored above the 

surface. The bottom exerts on the top a force —F^ , which 
by itself would create a displacement 



uf\x)= / Xi^x-x') M'V)) dx' (6) 
J —00 

This displacement is calculated in a way similar to the 
previous one (u*- 1 -*). As explained by Landau[^|, it involves 
a response function \ij which is the transposed of Xiji 
because the medium now sits above instead of below. One 
has to substract an energy 



60 



(2) _ 



1 



-j-00 



dxuf ] (a) l-F^ (x) 



(7) 



The cross terms cancel out between x and X an d we are 
left with a simple result for the elastic energy at this stage: 



[3 cl = I3 { 



/+00 
dx dx 1 ^ xu ( 3 
-00 x 



x-x')a { °> (x) a { °' (x) 



(o) 



(8) 

As explained in the appendix, we know X ij, so that we 
can solve this elastic problem analytically. Note that we 
could also use the method of images if we added coun- 
terforces 2^°' along the surface only. Comparing the two 
approaches is a useful check of a rather cumbersome alge- 
bra. 

The calculation is most easily carried out in Fourier 
space, according to the definitions 



fix) 



+ OO 



/ (k) e lkx dk 



1 r + °° 

f(k) = — f(x)e- ikx dx 

27T J-00 



(9) 
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For a periodic force Uj (k) = Xij (k) Fj (k). It is easily ver- 
ified that 



XH (x) = - 



The elastic energy is 

/3 cl = l3 {o) 



-Hex:- 



Xij (k) 



3 lkx dk 



(10) 



2tt 



dk 



(11) 



We only need carry integrals! 

But this is not the end of the story. We have not yet in- 
cluded elastic deformations due to the existence of the sur- 
face tension, that is to capillarity! The normal stress a zz 
at the surface is not zero. The surface is under mechanical 
equilibrium because the surface stress is compensated by 
the Laplace force <P Z which is exerted by the surface on 
the inside. It is a consequence of the curvature in the step 
profile: 

(12) 

That force creates an extra displacement it cap which in 
turn modifies the profile. We thus write in reciprocal space 

(k) = - 7 fc 2 «> + uT p ) 

u c z ap (k) = Xzz (k)<P z (k) 



whose solution is 



,cap ^ _ 



l + otk 

2, 



"fe = lk Xz Z (k) 

Because the displacement it cap modifies <P Z (k), the usual 
counterforce argument does not apply and it is simpler to 
write the net capillary energy as the sum of its bare value 
and the elastic energy stored in producing u cap 



□cap 



2n 



fo ° f ik 2 

dkl2_\ u f {k) 



cap 



(*0| 



.cap 



2 X zz (k) 



(13) 

(minimizing with respect to it cap reproduces the expres- 
sion of it cap ). We thus find 



cap 



2tt 



dk 



-fk 2 



2(1 + a*) 



(14) 



As expected, the elastic relaxation reduces the original 
capillary energy. 

Our program is now clear: first we calculate (k) 
and x (k), then we minimize (3 = /3 cl + /3 cap with respect 
to the depth h. 

Finally we calculate the step energy, which is the to- 
tal energy of the dislocation. The algebra is done in the 



Appendix. It leads to 

,2 




4ir (1 - <t 



— e 



-2kh 



1 + 2kh + 2k 2 h 2 - 2kd, 



{l + khf 
l + kd 



(15) 



where d a = (1 — cr)j/(i is our characteristic length. The 
term 1 is the logarithm in [3^ , the first three terms in the 
square bracket are the contribution of counterforces and 
the last one the effect of capillary forces. The denominator 
(1 + kd Q ) reflects the self consistent calculation of it cap . 
Note that capillarity is hidden in the single length d . We 
rewrite the curly bracket as 



1 



-2kh 



2khe- 2kh 



1 _do 
h 



1 + kh 



1 - d fh 
1 + kd n 



(16) 



It is easily verified that this quantity is stationary with 
respect to h when h — d a . That result is an exact version 
of our original estimate. Note that minimization is the 
same for every k: as a consequence it will persist for two 
steps or for a regular array. Altogether the total energy 
takes an incredibly simple form 







4tt(1 - a) 

,2 



1/T ° dk 



-2kd 



4tt(1 



In 



2d r 



(17) 
(18) 



where r w 0.577 is Euler's constant. All complications 
have disappeared! In view of the messy intermediate al- 
gebra there is probably a reason for that, but we do not 
know it. 

Eq. ( |l8| ) is our first conclusion and it is non trivial: 
it shows that the energy scale for steps is not 7 a, but 
/Lta 2 , down by a factor //a/7 ~ a/d . This will affect 
the roughening transition temperature of vicinal surfaces 
drastically. Step burying is physically important. 



3 The interaction between steps 

The interest of Fourier space calculations is the fact that 
they are immediately extendable to several steps. For a 
pair of steps a distance d apart, the strain and stress ten- 
sors are simply multiplied by a factor (l+e lkd ). The result- 
ing energy integral for two steps is modified by a factor: 



II 



ikd\2 



2 (1 + cos kd) 



(19) 



Note that the depth h = d is unchanged, as minimization 
was achieved for each k. The interaction between steps is 
contained in the term 2 cos kd. Using the above result for 
the energy of an isolated step we find that the energy for 
the two steps becomes: 



2f3(d) 



fia 



2tt(1 -a) 



1 /' r ° dk 

dk {l 



-2kd„ 



[1 + COS kd] 
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/id 



2tt(1 - a) 



r + ln 



2d, 



1, (■> 4d l 
2 ln 1 + # 



(20) 



The interaction energy is contained in the last term and 
for large d it is: 



£12 



(1 - a) d 2 



(21) 



It is interesting to compare this global repulsion with 
the successive results one would get at various steps of an 
analysis using image dislocations. 

(i) Using only image dislocation without compensating 
forces 2F( ' leads to the interaction energy of four dislo- 
cations arranged in a rectangle with height 2h and base d. 
We want half the total elastic energy e of these four dis- 
locations since we integrate the elastic energy over a half 
space. The starting point is the interaction energy of a 
pair of parallel edge dislocations, with Burgers vectors b\ 
and 62 pointing in the same direction, but either parallel 
or antiparallel. That interaction is given in the literature 

rani 



£12 



\ib x b 2 



lit (I -a) 



In 



R 



sin 



(22) 



where R and 9 are the polar coordinates of 2 with respect 
to 1 (polar axis along b) , and -R max an upper cutoff which 
will disappear for dislocation dipoles. The first term in the 
bracket yields a radial force which is repulsive for b\ = b 2 
but attractive if 61 = — b%. The dislocation 1 is repelled by 
its neighbour (dislocation 2), but it is attracted by the im- 
age 2' of this neighbour. We also recover the known result 
that a double dislocation with Burgers' vector 26 tends 
to split into two single dislocations b. The second term in 
the bracket of Eq. (g2|) yields a torque. For our rectangu- 
lar configuration, and burgers vectors of amplitude a, we 
find 



4 In- 



1 

2 4tt(1 



41n- 



Rn 



4 In- 



Rn 



4 In- 



Rr, 



2h 

. cos a 



cos 2 a 



In the bracket the first term is the energy of individual 
dislocations, the second the interaction of genuine dislo- 
cations 1 and 2 with their respective images 1' and 2', 
the third (1,2) and (l',2'), the last one (1,2') and (l',2). 
The angle a is the tilt of cross links, tan(a) = 2h/d. That 
result holds for any value of h/d as long as h and d are 
much larger than a. The elastic interaction between the 
steps corresponds to the last two terms in the bracket 



fj,a 



2tt(1 



[sin 2 a + ln(cos a)] 



(23) 



For small a, i.e. large distance d, the bracket is a 2 /2: the 
interaction is attractive, proportional to 1/d 2 as befits a 
step interaction. It turns into a repulsion for larger a, the 



change corresponding to the minimum of ei at a = tt/A, 
i.e. d — 2h. If true, such a behavior would lead to the 
prediction of a sharp angular matching at the edge of a 
fiat facet §. 

(ii) The next approximation is the inclusion of counter- 
forces at the surface. One finds that the 1/d 2 attraction 
is exactly cancelled out. The net step interaction is con- 
trolled by the next order in the expansion: it is a 1/d 4 
attraction. 

(Hi) But capillarity ruins the above result. As any relax- 
ation mechanism it induces a 1/d 2 repulsion that super- 
sedes the 1/d 4 term. We conclude that step interaction 
is controlled by the elastic strain field which is produced 
by the surface tension. Just as in the classic picture by 
Marchenko and Parshinjl(| it is repulsive and ~ 1/d 2 . 
Our picture does not depart from standard wisdom. In- 
terestingly, Eq. ( pl| ) can be written as 



£12 



(1 - ah 2 a 2 
n /id 2 



(24) 



which looks like Marchenko's result with a force doublet 



7a on the step; however, 
tributcd on a length do- 



the force doublet is /xa dis- 



4 The energy of a vicinal surface 

In practice we are not interested in a single pair of steps, 
but in a regular array with period d, describing a vicinal 
surface with average tilt <f> such that tan(/) = a/d. Due to 
that tilt the Burgers vector of dislocations is not exactly 
perpendicular to the real surface, which strictly speaking 
makes the elastic calculation far more complicated. We 
can ignore such a complication if we assume again that 
the lengths h and d are both much larger than the atomic 
spacing a: the angle (j> is then small and we retain only the 
lowest order in a/d and a/h while making no assumption 
on h/d. Once more, the depth h is the same as for a single 
step. 

We want to calculate the total energy per step ft and 
per unit area U = 8/d. We can proceed in the same way 
we did for two dislocations. Given the Fourier transform of 
the stress tensor of a single step er(°)(fc), we can generate 
the periodic superposition as follows: 



+ N 



i kdn Jo) 



(k), 



(25) 



where 2N + 1 is the number of steps. If the total energy 
of a single step is 



0i= b(k)dk 



(26) 



the total energy per step in a periodic arrangement be- 
comes: 



+N 

3 = V / dk 

' 2N+1 ^ J 



D i kd{n — 7n) 



b(k) 



n,-m— — N ' 
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2n 
~d 



+ 00 

E 



V— — 00 



2nv 



(27) 



Using the expression for b(k) given in Eq. ( fL7j ) we obtain: 
■2k,I.. (1 - e -^old) 



(3(d) = 



4tt(1 - a) 



d/(2irr 

E 



(28) 

The energy per unit area is U = (3(d) /d. The single term 
v = yields the capillary energy 7a 2 /2d 2 associated to 
the surface lenghtening due to an average tilt (f>. 

The summation can be done analytically for d /r 
large and it gives: 



(3(d) 



4tt(1 - a) 



2ird 



T + ln 



27rr Q 



In 1 



3 -4ird /d 



for 7, it is the value of the surface stiffness at Tr. In the 
case of interest here, where the lattice periodic potential 
is a small perturbation ("weak coupling"), one can neglect 
the critical variation of 7 near Tr and use a temperature 
independent value for 7. 

We are concerned here with facets of a vicinal surface 
(1, 0,p), which correspond to an arrangement of "primary" 
(1,0,0) steps commensurate with the in plane periodic- 
ity of the underlying crystal. (If the crystal was periodic 
in one direction only, it would be a stack of structurless 
plates, and there would be no (1,0, p) facet.) This prob- 
lem is discussed in the accompanying paper 17fl. The use 
of the universal formula ( |3l| ) needs a little care: 

(i) a becomes a(p) — a/p. Indeed, for a cubic crystal, 
one goes from one (1,0, p) plane to the next by translating 
primary steps by a lattice spacing a. 

(ii) 7 is anisotropic and it should be replaced by the "sur- 
face stiffness" 

<9 2 7 

Tij = 7% + an an ( 32 ) 



dOidO; 



where we used YZ=i l l v = r + lni/ m + 0(1/ v m ) and After diagonalization of 7^, Eq. (|31|) is replaced by 
Y^ =l e~ vX jv = — ln(l — e~ A ). We rearrange the expres- 



sion in order to separate the core dependence. We finally 
find: 



kBT R (p) = - y/jxly a(pf 

7T 



(33) 



(3(d) 



jia 



An(l - a 
2d, 
r, 



r + In — 



2itd n 



In 



I _ e -iivd /d 

Am d Q I d 



.(29) 



As mentioned above the relevant scale for the d depen- 
dence of (3 is d Q . The expansion of (3(d) for large d begins 
with a positive 1/d 2 term, indicating repulsion of steps: 



(3(d) = 



4tt(1 - a) 



r + ln(^ 



2vr 2 d ; 



2^2 1 



3d 2 



(30) 



The step-step interaction coefficient is correctly z~2^=i l/^ 2 
= 7r 2 /6 times the two-step interaction energy given in (f2l|). 
For small d the leading term is simply the capillary energy 
needed to lengthen the surface in order to tilt it ~ 1/d. 



5 Roughening temperatures 

It is now well known[|||l^] that, as their temperature de- 
creases, crystals are facetted in an increasing number of 
different directions. In any rational direction, there is a 
roughening transition at a temperature Tr. Below Tr, the 
surface is "smooth", i.e. facetted or flat. It is localized by 
the influence of the crystal lattice. Above Tr, the surface 
is "rough" like a liquid surface and rounded on the equi- 
librium shape of the crystal. The roughening temperature 
is given by the universal relation |15[ |l2|| 



The quantity j x is controlled by step compression and it 
is given by: 

1 d 2 U 

^ = TTT ( 34 ) 
a z dn 

As for j y , it appears in the energy cost for an undulation 
of steps which move in the x direction by an amount u = 
uoexp(ik y y). Whenever a step sweeps by a given point 
the interface is raised by a, from which we infer a height 
change Sz = na a u. If each step rotates by an angle Q = 
du/dy, there is both lengthening and compression of the 
steps. Since one has to replace the energy element (3(n)dx 
by (3(nj cosQ .(dx / cosQ , the energy change is of order £ 2 , 
and writes 



C 2 

5U =-2 



2 d(3 



n(3 + n — 
dn 



from which we infer 



7y 



- 1 fc 2^ 2 

2 dn ~ 2 lyk y dZ 



1 dU 



dn 



Alltogether, we find 



U'U" 



(35) 



(36) 



(37) 



kBT R = -70 

7T 



(31) 



where a the periodicity of the surface energy along the 
normal to the interface, in practice the step height. As 



Since the interaction at small angle, i.e. small n, is pro- 
portional to 1/d 2 , the asymptotic form of U(n — > 0) is 
n(3 + 5n 3 and the effective stiffness writes \/6(3d /a 2 in this 
limit. It is independent of the tilt angle or step density. It 
only involves the product of the step density by the step- 
step interaction coefficient S that can be inferred from 
©■ 

We are now equipped to infer Tr from the preceding 
section. Note that we do not need to expand in powers of 
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n, we can calculate Tr for arbitrary values of n do, the only 
assumption being d > o as before. The d dependence of 
the factor yJU'U" /n is not strong and can be conveniently 
studied by defining the dimcnsionless Tr: 



Tr = -^a(p) 2 T R . 



(38) 



scale of Tr. As a consequence of the small value of the 
elastic modulus /i, the step-step interaction is large, but 
in the expression for the roughening temperature, it ap- 
pears multiplied by the step energy which is small. We 
have shown that the product of these two quantities is 
nearly independent of ji. 



The d — > oo limit of Tr involves only the step energy and 
the step-step interaction: 



T R (d 



°° ) = 1 2\r( 2 -t 



r 



(39) 



It thus depends on the core radius r Q . The next order in 
the development in 1/d is quadratic and it can change 
sign as a function of d /r due to the fact that the 1/d 4 
term in the interaction energy of two steps is negative [cf. 
Eq. @]. In the limit where d -> 0, T R = 1. The full 
expression coming from (^9|) interpolates between these 
two limits and it is shown in Figure |l| for different values 
of the ratio r /d . It indicates the departure from a simple 
quadratic variation of Tr with the step density n = 1/d. 



1.2 



|: 



0.2 



0.4 

d 



, 0.6 

d 



O.J 



Fig. 1. Tr for different values of the parameter d /r as 
a function of n = d /d. The parameter d /r is 4, 10, 100. 
starting from the bottom. The dashed line corresponds to a 
d /r = 4 which is a reasonable value for the experiment of 
Ref. §. 



The crucial result is that Tr is of the order 7 a 2 , de- 
spite the fact that elasticity played a central role. The 
lower energy scale /ia 3 is compensated by the lower density 
scale l/d : since the effective stiffness has dimension U/n, 
elasticity drops out. That might be seen as trivial! But 
a naive argument bared in the usual TLK model would 
predict a step energy (3 ~ 7 a and an elastic interaction 
between steps 8 ~ (7a) 2 /fid 2 . The resulting effective 7 
is then 7(7/^ a) 1 / 2 : an exact calculation does reduce the 



6 Comparison with the experimental shape of 
lyotropic crystals 

Let us finally consider the number of facets observed by 
Pieranski et al.jl],p|]. These crystals show 6 beautiful stair- 
cases of facets around large (1,1,2) facets which we con- 
sider as the main ones (see fig. 5 in ref.^J). These crystals 
are cubic with Ia3d symmetry, which explains why the 
main facets are not (1,0,0). The last facet, i.e. the one with 
the highest Miller indices, is (9,8,15) (see fig. 7 in ref.Q) at 
293K, the temperature of the experiment. It corresponds 
to the smallest tilt angle with respect to the (1,1,2) facet, 
namely 4 degrees. Pieranski et al. explain that the distance 
between equivalent planes in a [h,k,l] direction is given by 



dhkl — 
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y/h 2 + k 2 + P 



A 



(40) 



This means that the "primary" step height is a — 48. 2 A 
on (1,1,2) facets. As for the "secondary" steps on vici- 
nal facets, one finds respective heights of 7.02, 6.13, 5.44 
and 4.89 A for (8,7,13), (9,8,15), (10,9,17) and (11,10,19) 
facets. 

The large value of a is responsible for the large number 
of facets. Indeed, according to Eq. (pl|), the roughening 
temperature of (1,1,2) facets is very high, 26800 K. We 
now understand that room temperature (293K) is very 
cold compared to the energy scale involved in the rough- 
ening of the main (1,1,2) facets. If we ignored the correc- 
tion Tr in Eq. (|38|), we would predict that the last facet 
at 293K corresponds to a step height slightly larger than 
a/^26800/293 = 5.04A. This would allow (10,9,17) facets 
to exist at 293 K. 

In the case of their crystals, we find a length do = 
3.46a = 167A. According to our calculation (see figjl]) 
the correction Tr is about 0.84 for tilt angles of order 4 
degrees, i.e. do/d = 0.24, so that we predict that the last 
facet should have a step height larger than a minimum 
value 5.50 A. This is just enough to avoid the existence 
of (10,9,17) facets, and allow that of (9,8,15) facets, as 
seen by Pieranski et al. Equivalently, we could say that 
our predictions for the roughening temperatures are 364K 
for (9,8,15) surfaces which have thus to be facetted at 
293K, but 287K for (10,9,17) surfaces which should thus 
be rough at 293K. 

Given the number of approximations in our theory, the 
agreement with the experiment is surprisingly good. The 
calculation is made for a simple cubic lattice where the 
Burgers vectors are normal to the main (001) faces. Since 
Pieranski et al. studied lyotropic crystals which have a 
laid symmetry and tilted Burgers vectors, our model is 
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clearly oversimplified. We believe that two things should 
remain if one introduced such details in a more accurate 
calculation: 

1- The highest order facet has to correspond to a step 
height of order one tenth of the main one because Tro is 
about 100 times 293 K. 

2- The correction Tr being smaller than one, one observes 
a slightly smaller number of facets than what a simple 
quadratic law would predict. Eventually, it is interesting 
to explain why we did not consider entropic interactions 
between steps. The origin of this other type of interaction 
is the no-crossing condition of steps |H| . Its amplitude has 
been calculated fl4| and measured [p~3|T It is 



6 



{k B Tf 1 
f3 d 2 



(41) 



Given the low value of the step energy that we have found, 
one could think that the entropic interaction is a signif- 
icant effect. In the case of soft crystals, the ratio of the 
entropic interaction to the elastic interaction is indepen- 
dent of the shear modulus fi. It would be large near the 
first roughening temperature, but it is very small at tem- 
peratures two orders of magnitude lower. 

We are very grateful to Pawel Pieranski and Paul Sotta, 
who introduced us to the physics of their crystals and com- 
municated their last results before publication. 



A Solution of the 2D elastic problem 

In this appendix we provide details on the solution of the 
two dimensional elastic problem (every quantity depends 
only on the depth z and the coordinate x). We consider 
a half space isotropic elastic medium z < 0. The basic 
ingredient is the response function relating the displace- 
ment Ui (x) at the free surface to the force Fj (x r ) applied 
elsewhere on the same surface 



Ui (x) 



Xij (x - x') Fj (x 1 ) dx' 



(42) 



Xz 



1 — a 



Xxz 



1 - 2a 
2iku 



(47) 



That same result could be obtained from the response to a 
point force given in Ref . || , integrating first on the dummy 
coordinate y and then Fourier transforming. 

The displacement and stress fields due to an edge dis- 
location along the y axis in an infinite medium are well 
known 



» - 



,C°) 



b 

2^ 



b 
2^ 



x 

arctan — -+ 

z 

1 -2a 



1 



4(i-*; 

lib x 



2(1-0") x 2 + z 2 
log (x 2 + z' 



1 



2{l-a)x 2 + z 2 



(32 



2vr(l-a) (s2 +s 2)< 
lib z(z 2 -x 2 ) 



2w(l-a) ( X 2 +J 2) 2 

~z is the distance from the core. If the dislocation is buried 
at a depth h below the origin, z = z+h. The corresponding 
stress and displacement at z = are in Fourier space 



(°) (]?) - ^ b 
zx [ >~ An(l-a) 



k\ h exp [— |fc| h] 



A") 



(*) 



^tegn(k) + exp { _^ h] 



4tt(1 - a) 



W - "BE 



\k\h 



2(1 



exp [— \k\ h] 



From that we infer the extra displacement iti (fc) due 
to relaxation of counterforces, in which we identify the 
contributions due to F z °^ and to fJ°^ 



(k) = - X zx (k)<>(k) 



(°) (h\ - 



6(1 - 2a) 
8irik(l - a) 



\k\he' Wh 



{k) = - Xzz (k)a^ ( k ) = - — [l + \k\h] 



,-\k\h 



The problem is solved via Airy functions (see ref.§) The net surface profile as controlled by elasticity is 



8 2 A 



dz< 



d 2 A 
dx 2 



d 2 A 
dxdz 



(43) 



(48) 



For a single Fourier component the biharmonic function 

A (x, z) is 



It is easily verified that uf = 2-iti 1 ^ : the shape of the 
surface is controlled by by F^ only. Such a very simple 
result is not an accident: instead of using a single disloca- 
tion with two counterforces Fz^ and Fx , we could have 
used an image picture with a doubled tangential coun- 
terforce 2Fx ■ The displacement would then read uf = 

2uf^ + 2u z lxS> ; identification of the two expressions of uf 
yields the above result! (As mentioned in the text explicit 
comparison of the approaches at every stage of the calcu- 
lation is an useful trick for hunting sign mistakes). 

In addition to the surface profile, we need the net elas- 
tic energy which we write as 



A (x, z) = (a + bz) exp z + ikx] 



(44) 



The constants a and b are fixed by the surface conditions 
a zl (0) = Fj, yielding 



-F z /k 2 



b = iF x /k + F z /\k\ 



(45) 



From <7jj we infer the strain then the displacement 
u x integrating u xx . We thus know du x /dz, from which we 
infer du z /dx : u z follow upon integration over x. The final 
result in momentum space is 



Ui (k) = Xij (k) Fj (k) 



(46) U cl = U { 
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4tt(1-(7) Jo 



Xzz (k)a^(k) +Xxx(k) 
1/r ° dk 



(k) 



k 



{1- [1 + 2kh + 2k 2 h 2 ] e~ 2kh } 



The infrared divergence at k — cancels out (the upper 
distance cut off is h instead of -R ma x), but the large k diver- 
gence remains: it describes the dislocation core! The last 
step is inclusion of the capillary correction, as described 
in the main text. The relevant quantity is 



ctk = yk 2 Xzz (k) = kd 



(49) 



Since the relevant k are w l/d , the need for a selfconsis- 
tent calculation is apparent. The net capillary energy, as 
"screened" by elastic response, is 



[/cap 



= 2tt 




uf(k)+uT P (k)\ + 



uf(k)\ 



.cap 



(k)Y 



2Xzz (k) 



The expression for the total energy used in the text follows 
at once. 
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